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I. 


Introduction. 


N-isopropylcarbazole  (NIPC)  is  a  pyroelectric  molecular  crystal 
1  2  3 

exhibiting  strong  tribo-  and  pyroluminescent  properties  associated  with 


the  polar  structure  of  the  crystal.  A  structural  phase  transition  in  NIPC  at 

,21 


ca.  137  K,  taking  place  between  the  C^v  (iba 2)  high  temperature  structure 

C 

and  the  (Pbc21 )  low  temperature  structure,  has  been  reported  in  the 


literature4*^.  Dilatometric,4  pyroelectric,2,4  and  specific  heat6  measure¬ 


ments  reveal  the  first  order  character  of  this  transition.  The  enthalpy  of  the 


-1 


transition  is  found  to  be  AH  =125  cal  mol  ,  a  value  higher  than  most  dis- 


7  8 

placive'  or  even  order-disorder  phase  transitions  in  molecular  crystals.  Cal¬ 


orimetric  measurements0  reveal  an  additional  anomaly  of  the  specific  heat  around 
185  K,  indicating  the  possible  existence  of  another  (non-first  order)  phase 
transition. 

Very  little  is  known  about  the  mechanism  of  the  NIPC  phase  transition 
at  137  K.  From  the  symmetry  point  of  view,  NIPC  is  a  rare  example  of  a 
molecular  crystal  exhibiting  a  phase  transition  without  any  point  symmetry 
breaking;  only  the  translational  symmetry  changes  at  the  transition  (the 
volume  of  the  primitive  unit  cell  doubles  in  the  low  temperature  phase). 


Chloranil^*^  is  another  well  known  example  of  such  a  nonferroic^1  phase 


transition  among  organic  solids. 


A  group  theoretical  analysis  of  the  transition  in  NIPC  suggests  that 

.5 


the  transition  takes  place  at  the  Brillouin  zone  boundary  and  that  the  order 
parameter  transforms  according  to  a  one  dimensional  representation  of  the 
group  of  the  wavevector  (only  one  vector  of  the  star  of  the  wavevector  is 
involved  in  the  transition).  Bilinear  coupling  of  the  order  parameter  to 
macroscopic  quantities  at  the  zone  center  (e.g.,  strains)  is  thus  forbidden 
by  symmetry.  The  lowest  order  coupling  permitted  by  symmetry  is  then  a  third 
order  nonlinear  interaction.  Although  this  kind  of  transition  does  not 


involve  new  macroscopic  tensor  components  below  the  phase  transition, 
anomalies  of  the  elastic  constants  may  arise  due  to  the  higher  order 
couplings,  as  is  the  case  for  chloranil. ^ 

The  molecular  mechanism  of  the  phase  transition  in  NIPC  is  not 

5 

completely  clear.  Recent  crystallographic  data  have  shown  large  librational 
motion  of  the  NIPC  molecules  at  room  temperature.  In  particular,  libration 
around  the  molecular  axis  of  largest  moment  of  inertia  exhibits  surprisingly 
large  amplitude  (7.1°),  at  least  twice  the  amplitude  of  the  two  other 
principal  axis  librations.  The  amplitudes  of  the  librations  become  almost 
equal  and  significantly  smaller  in  the  low  temperature  phase  and  only  small 
rotation-translation  shifts  of  the  molecules  occur  at  the  transition.  The 
carbazol  group  of  the  molecule  itself  is  planar  in  both  phases.  The  only 
significant  change  in  molecular  geometry  at  the  phase  transition  is  a  rotation 
of  the  isopropyl  group  around  the  C-N  bond  by  an  angle  of  14°  and  10°  for 
layers  I  and  II,  respectively.  These  results  show  that  the  phase  transition 
in  NIPC  does  not  affect  the  molecular  conformation  in  an  essential  way:  the 
phase  transition  may  be  considered  as  an  order-disorder  or  rotation-translation 
coupling  type  transition.  Rotation-translation  coupling,  as  well  as  "freezing" 
of  the  orientational  motion  of  molecules  in  the  low  temperature  phase,  has  been 
previously  reported  to  characterize  phase  transitions  in  such  molecular  crystals 
as  s-triazine,^  benzil,1^  and  chloranil.^® 

As  was  already  mentioned,  the  NIPC  crystal  is  pyroelectric,  ar.d  thus 
also  piezoelectric,  in  both  phases.  The  molecular  dipole  moments  are  directed 
along  the  C-N  bonds  of  the  molecules  which  are  oriented  about  the  direction  of 
the  two-fold  c-axis;  the  spontaneous  polarization  of  the  crystal  is,  therefore, 
parallel  to  this  axis.  Due  to  the  piezoelectric  coupling,  dielectric  and 
elastic  anomalies  cannot  be  separated,  in  general.  Therefore,  some  of  the 
observed  elastic  constants  must  be  corrected  by^ 
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ca6  *  cag  "  ekag/Ek 


in  which  cag  is  a  bare  elastic  constant,  0  is  a  piezoelectric  constant, 


kng 


and  is  a  bare  dielectric  constant.  Since  the  anomalous  part  of  the 

4 

dielectric  constant  of  NIPC  is  much  smaller  than  its  normal  part,  the  main 
effect  of  the  piezoelectric  coupling  is  to  shift  the  absolute  values  of  the 
elastic  constants.  Ignoring  this  effect,  interactions  between  the  order 
parameter  and  strains  can  be  separated  from  interactions  between  the  order 
parameter  and  polarizations. 

In  this  report  we  present  a  complete  study  of  Brillouin  scattering  in 
NIPC  at  one  atmosphere  pressure  between  295  and  ca.  100  K.  The  elastic  constant 
tensor  is  determined  and  the  temperature  dependences  of  the  elastic  constants 
are  given.  The  observed  anomaly  of  the  a-axis  longitudinal  acoustic  phonon 
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associated  with  the  e1  strain  is  discussed  on  the  basis  of  Landau  mean  field 
theory  and  interpreted  in  terms  of  a  Debye  relaxation-dispersion  of  the  elastic 
constant. 

II.  Experimental  Procedures. 

Single,  transparent  crystal  boules  of  NIPC  (ca.  5  cm^  in  volume),  exhibit¬ 
ing  a  [lOO]  perfect  cleavage  plane  and  [010]  non-perfect  one,  are  grown  from  the 

1 5 

melt  using  a  Bridgman  technique.  Oriented  crystals  are  cleaved  with  a  razor 


| 


blade  and/or  cut  with  a  wire  saw  into  parallelepiped  shaped  samples  (ca.  4  x  4  x 
3  mm).  The  saw-cut  faces  are  polished  on  a  tissue  soaked  with  ethyl  acetate. 

In  order  to  reduce  surface  scattering  and  prevent  sample  sublimation  during 
experiments,  crystal  faces  are  coated  with  a  thin  layer  of  silicon  oil  prior  to 
the  experiment. 

Fight  angle  scattering  experiments  are  performed  using  three  different 
kinds  of  samples  characterized  by  different  orientations:  1)  faces  perpen¬ 
dicular  to  the  three  crystallographic  axis  -  this  orientation  allows  for  measure- 


ments  of  phonons  in  the  directions  [lio],  [ 1 01 ] ,  and  [oi 1 ] ;  2)  faces  oriented 

45°  with  respect  to  the  crystallographic  axes  in  the  ab  plane  -  phonons  in  the 

directions  [lOO]  and  [ 01 0 ]  can  be  measured  in  this  orientation;  and  3)  faces 

oriented  45°  with  respect  t.o  the  crystallographic  axes  in  the  ac  plane  -  either 

[lOO]  or  [ 001 ]  phonons  can  be  obtained. 

The  Brillouin  scattering  apparatus,  including  measuring  techniques  end 

12  13 

temperature  control,  used  in  these  experiments  has  been  described  previously.  ’ 

A  Spectra  Physics  model  165  argon  ion  laser,  operating  at  5145  X  (200  mW)  in  a 
single  cavity  mode  is  used  for  the  scattering  source.  Spectra  are  obtained  with 
a  Burleigh  Fabry-Perot  interferometer  system  in  triple-pass  operation  with  a 
standard  finesse  of  45  to  50.  To  analyze  the  spectra  and  to  record  the  tempera¬ 
ture  dependences  of  the  acoustic  phonons,  free  spectral  ranges  between  15  GHz  to 

30  GHz  are  chosen.  Typically,  50  or  100  accumulations  of  a  Is  duration  scan  are 
accumulated  and  averaged  in  a  multichannel  analyzer  and  computer. 

III.  Results. 

Employing  the  standard  procedure  for  Brillouin  scattering,  the  velocity 

1 5 

of  sound  v  in  an  optically  anisotropic  crystal  can  be  calculated  by 


2 it  /  2  2  _  \J< 
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in  which  n  ,n  are  the  refractive  indices  for  incident  and  scattered  light, 

X  8 

respectively,  X  is  the  incident  light  wavelength  in  vacuum  and  6  is  the 
scattering  angle  (90°  for  right  angle  scattering).  Thus,  knowing  the 
frequency  shift  (u»g)  and  the  polarizations  of  the  incident  and  scattered 
light  from  the  experiment,  one  can  easily  calculate  the  sound  velocity, 
providing  refractive  indices  of  the  crystal  are  known. 

The  refractive  indices  of  NIPC  are  measured  at  5145  X  by  a  prism 
minimum  deviation  method.  The  accuracy  of  this  method  in  our  case  is  only 


■jrrjrr sn 


slightly  better  than  5%}  the  error  is  mostly  due  to  the  poorly  determined 
geometry  of  the  prisms  and  the  quality  of  their  surfaces.  The  values  of  the 
refractive  indices  are  subsequently  verified  and  adjusted  by  comparing  their 
ratios  to  the  ratios  of  the  Brillouin  shifts  measured  along  the  principal  axes 
in  backscattering  geometry  for  different  polarizations.  Refractive  indices 


obtained  in  this  way  and  used  for  calculations  of  sound  velocities  are  n  = 

a 


1.846,  n^  *  1.648,  nc  *  1.664.  Elastic  constants  of  NIPC  are  calculated 

3  c 

using  the  density  P  =  1.15  x  10  kgm  taken  from  x-ray  data. 


A.  Determination  of  the  elastic  constants  at  room  temperature. 

Elastic  constant  tensor  and  piezoelectric  constant  tensor  for  both 


orthorhombic 

phases 

of 

NIPC  are 

given  by 
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respectively.  The  equation  of  motion  for  the  piezoelectric  crystal  is  of  the 
form 


<cijki  -  ov  stk>  \ 


(2) 


in  which 
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are  the  elastic  constant  coefficients  modified  by  the  piezoelectric  coupling, 
and  the  m's  are  the  direction  cosines.  The  solutions  of  the  equations  of 
motion  corresponding  to  the  phonons  observed  in  our  experiments  are  shown  in 
Table  I.  Additional  terms  of  the  form  of  eq  (3)  in  the  Table  describe  the 
change  of  ine  velocity  of  those  phonons  which  couple  to  the  electric  field  via 
piezoelectric  coupling.  As  already  mentioned  in  the  Introduction,  the  dielec¬ 
tric  properties  of  NIPC  do  not  show  a  large  anomaly  at  the  phase  transition.^ 

As  a  result,  interactions  between  the  order  parameter  and  strains  can  be,  in 
principle,  separated  from  interactions  between  the  order  parameter  and  polariza¬ 
tions.  Thus,  only  the  shifts  of  absolute  values  of  the  elastic  constants  will 
be  considered  in  the  determination  of  the  elastic  constant  tensor  and  its  tem¬ 
perature  dependence.  Values  of  these  shifts  cannot,  however,  be  calculated 
from  eq  (3)  due  to  lack  of  experimental  data. 

On  the  other  hand,  the  change  of  the  acoustic  mode  velocity  due  to  the 

17 

piezoelectric  coupling  can  be  estimated  from  the  electromechanical  coupling 

©ff*  2  IS 

coefficients  k.  specifically,  c  -  c/(l-k  ).  Preliminary  measurements  7 
*  J 

have  shown  that  k^  *  0.15,  k^  *  0.18,  k2^  *  0.11:  k ^  and  could  not 
be  measured  due  to  weak  piezoelectric  signals.  Employing  these  values  of  the 
electromechanical  coupling  coefficients,  the  largest  change  in  the  sound  veloc¬ 
ity  caused  by  piezoelectric  coupling  is  estimated  to  be  1.6/6.  Thus,  the  exper¬ 
imentally  observed  elastic  constants  c  differ  from  bare  elastic  constants  by 
less  than  -3$.  The  piezoelectric  coupling  will  not  be  included  in  the  ensuing 
discussion,  as  it  is  of  little  importance  to  the  elastic  anomalies. 

The  values  of  the  elastic  constants  of  NIPC  are  calculated  using  the 
appropriate  relations  from  Table  I.  c^  ,  c22,  C55*  and  c^g  are  obtained 
directly  from  X^ ,  X^,  X^,  X^  and  Xy,  X^  and  X^,  respectively,  with  an  accuracy 
of  ca.  2/6.  As  we  are  not  able  to  detect  the  transverse  mode  governed  by  the  c^4 


elastic  constant  directly,  its  value  is  determined  from  the  velocity  of  the  TA 
phonon  propagating  in  the  [lio]  direction  using  relation  X,^.  The  off-diagonal 
elements  of  the  elastic  tensor  c^*  and  c^^  are  calculated  from  the  rela¬ 
tions  Xg  -  Xg,  X1 1  -  X12  and  X^  -  X^,  respectively.  Their  values  are 
determined  with  a  relatively  large  error  i'30%)  due  to  the  combination  of  the 
elastic  constants  in  the  employed  equations,  as  well  as  to  the  relatively  small 
values  of  the  off-diagonal  elements.  The  room  temperature  set  of  elastic  con¬ 
stants  found  in  the  above  manner  is  given  in  Table  II. 

Using  the  values  of  the  elastic  constants  from  Table  II,  angular  sound 
velocity  diagrams  can  be  calculated.  These  diagrams,  given  in  fig.  1,  reflect 
the  anisotropic  nature  of  the  crystal  structure.  The  c^  elastic  constant 
which  governs  the  velocity  of  the  [lOO]  LA  phonon  is  relatively  large,  whereas 
the  molecular  interactions  in  this  direction,  which  is  perpendicular  to  the 
cleavage  plane,  are  expected  to  be  weak. 

B.  Temperature  Dependences  of  the  Elastic  Constants. 

The  temperature  dependences  of  the  six  diagonal  elastic  constants 
are  shown  in  figs.  2-4.  The  c ^ »  c44*  8nc*  c 55  elastic  constants  exhibit 
the  step-like  changes  expected  for  a  first  order  phase  transition.  Only  a 
small  effect  is  seen  for  the  temperature  dependence  of  c^*  An  interesting 
anomaly  around  the  phase  transition  at  137  K  is  found  for  the  c^  elastic 
constant,  which  involves  the  e^  strain  in  the  crystal  and  governs  the  propaga¬ 
tion  of  the  LA-phonon  is  the  [lOO]  direction.  Fig.  5  shows  the  temperature 
development  of  the  [lOO]  LA-mode  spectrum;  the  other  mode  seen  in  this  figure 
is  the  [lOO]  TA-mode  governed  by  the  c^  elastic  constant.  The  Brillouin  peak 
corresponding  to  the  [lOO]  LA-phonon  in  the  spectrum  exhibits  considerable 
broadening  at  the  phase  transition.  The  latter  phenomenon  typically  indicates 
a  higher  sound  attenuation  at  the  phase  transition,  possibly  due  to  the  critical 
fluctuations.  The  temperature  dependences  of  the  halfwidth  of  this  Brillouin 
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peak  and  its  intensity  are  shown  in  fig.  6. 

The  anonalous  temperature  dependence  of  the  c^  elastic  constant 
shown  in  fig.  2  may  be  fit  using  a  phenomenological  power  law  equation  of  the 
form 

C11(T)  =  C°1  -  aT  -  be~P  (4) 

in  which  the  linear  temperature  term  accounts  for  the  behavior  of  the  elastic 

constant  in  the  absence  of  a  phase  transition  and  the  last  term  represents  the 

form  of  the  critical  anomaly.  e  in  eq.  (4)  is  equal  to  (T-T  /T  )  or  (T  -T/T  ) 

for  the  high  and  low  temperature  phases,  respectively.  The  values  of  the  para- 

o  9-2 

meters  which  give  the  solid  lines  in  fig.  2  are  ^  *  15-11  x  10  II m  ;  a  =  1.5  x 
-2  -1 

10  K  ,  b  *  0.35,  Tq  *  135*5  K  and  p  *  0.67  for  the  high  temperature  phase,  and 

c°  *  17.19  x  109  Nm"2,  a  -  1.5  x  10'2  K"1  ,  b  -  2.05,  T  *  140.8  K  and 
1  1  0 

p  *  0.26  for  the  low  temperature  phase. 

Temperature  dependences  of  the  c^,  c^  and  elastic  constants  can¬ 
not  be  obtained  directly;  however,  they  may  be  estimated  from  the  temperature 
dependences  of  the  modes  propagating  along  the  diagonal  directions  [llOj,  [lOl] 
and  [oil]  (see  Table  II).  The  temperature  behavior  of  these  modes  is  given  in 
fig.  7.  Calculating  the  temperature  dependences  of  c^,  c^  and  C2^  accord¬ 
ing  to  the  previously  described  procedure  employed  in  the  determination  of  their 
room  temperature  values,  one  finds  that  the  off-diagonal  elastic  constants  do 
not  shown  any  particular  anomalous  behavior  at  the  phase  transition.  In  other 
words,  the  temperature  dependences  of  the  diagonal  elastic  constants,  c^  in 
particular,  fully  account  for  the  anomalies  shown  in  fig.  7.  The  temperature 
dependences  of  the  and  consequently  c^  elastic  constants  cannot  be 
determined  in  the  low  temperature  phase  due  to  cracking  of  the  sample  at  the 


phase  transition.  In  fact,  several  samples  are  used  to  obtain  most  of  the 


temperature  dependences  shown  in  this  paper. 

Elastic  constant  anomalies  which  could  confirm  the  existence  of 
another  phase  transition  at  ca.  185  K,  suggested  on  the  basis  of  calorimetric 
measurements,^  are  not  observed.  Nonetheless,  a  change  in  slope  of  the 
temperature  dependence  of  c^  is  found  near  185  K. 

IV.  Discussion. 

The  present  Brillouin  scattering  study  reveals  the  softening  of  the 
[lOO]  LA-mode  governed  by  the  c^  elastic  constant,  both  above  and  below  the 
phase  transition  at  137  K.  Simulataneously  a  -200$  increase  of  the  attenuation 
(halfwidth)  of  this  phonon  occurs  at  the  transition.  To  explain  these  experi¬ 
mental  findings,  the  Landau  mean  field  theory  and  dispersion-relaxation  proc¬ 
esses  will  be  considered. 

A.  Landau  free  energy. 
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Consider  the  group- theoretical  analysis  of  the  C^v  phase  transi¬ 

tion  in  NIPC.  In  the  Brillouin  zone  of  the  high  temperature  structure  (T  v), 

five  symmetry  lines  possess  the  Cgv  point  symmetry,  appropriate  to  the  consid- 

Id  2 

ered  transition;  however,  only  three  points  of  this  symmetry  k  * 

— E  a 

0,0),  ■  (0,  g — ,  0),  and  *  (0,0,  - — ),  all  located  at  the  Brillouin 

zone  boundary,  give  doubling  of  the  primitive  unit  cell  volume  at  the  transi¬ 
tion.  Thus,  the  transition  takes  place  at  the  Brillouin  zone  boundary.  As  the 
star  of  the  wavevector  has  only  one  arm  and  is  non-degenerate,  the  order  para¬ 
meter  (e.g.,  optical  phonon)  driving  the  transition  belongs  to  a  one  dimensional 
irreducible  representation  at  the  zone  boundary. 

1 9 

Employing  standard  group  theoretical  methods  ,  the  Landau  free  energy 
may  be  expressed  in  terms  of  the  order  parameter  as 


9 


G^  =  ','?A(T)n2  +  1/4Bn4  +  1/6Dn6  ♦  ... 


(5) 


in  which  A(T)  =  a(T-To)  stands  for  harmonic  restoring  forces  of  the  optical 
mode.  The  elastic  part  of  the  energy  is 


G  ,  =  1/2  1.  .c.  .e . e . 
el  '  ij  lj  i  j 


(6) 
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in  which  only  the  term  Gg^  =  1/2c^e^  will  be  considered  because  only  the 
c^  elastic  constant  exhibits  anomalous  behavior  in  NIPC. 

As  the  bilinear  coupling  between  the  zone  boundary  order  parameter 
and  the  zone  center  strains  is  symmetry  forbidden,  the  lowest  order  couplings 
which  may  be  considered  are  a  third  order  nonlinear  interaction  and  a  biquad¬ 
ratic  coupling.  Thus,  the  free  energy  interaction  term  is  represented  by 


=  -1/2(Ke.n2  +  Le2n2)- 
e  n  1  i 


(7) 


The  total  Landau  free  energy  is  then  of  the  form 

G  =  1/2A(T)n2  +  1/4Bn4  +  1/6Dn6  ♦  1/2cne2  -  1/21^  p2  -  1/2Le2p2.  (8) 

The  solution  of  the  Hamiltonian  including  only  the  third  order  nonlinear 

coupling  is  already  known  in  the  literature  for  the  phase  transitions  in  terbium 

molybdate2^*  and  chloranil^®.  The  theory  predicts  a  continuous  downward  displace 

ment  of  the  elastic  constant  below  the  transition,  whereas  above  the  transition 
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the  elastic  constant  is  unaffected  by  the  coupling.  In  the  case  of  NIPC  the 


solution  has  the  form 


■xrziTFjr^ir^Hxs^ir*.  irw  ir*  tf^w *  \r*  ir"  ir*  \r*  r*\r*  \r*\r*  g^gw  py  VTnrwB*  v*v*v*  it*  y^yt  y?  yg  yy  t»  v  *y  yw  v- 


/c  *  (T>V 
1  11 
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:,  ,  -  - - - - -  ,  (T  <  T  ) 

11  2B  +  4Dn2 


(9) 


in  which  fj  represents  the  spontaneous  value  of  the  order  parameter  in  the  low- 
temperature  phase. 

Including  both  linear-quadratic  and  biquadratic  interactions,  as  in 
eq.  (7),  the  equilibrium  values  of  the  order  parameter  (n)  and  strain  (6^  are 
obtained  by  setting  the  generalized  forces  equal  to  zero 


!i  ■  (%)  -  m ° 


(10) 


yielding, 


A(T)  +  (  B  - 


o  -2 

'n  - L" 


)  n2  +  (  D  - 


LK 


.4 


)  n  -  0  (11) 


,  0  t'2n 

(cn  -  Ln  ) 


and  (cn  "  Lp^)  “  l/2Kn^  =  0. 


(12) 


Eq.  (12)  has  been  used  to  eliminate  e1  from  eq.  (11).  In  order  to  ensure  that 

2 

the  transition  is  of  the  first  order  (which  is  appropriate  for  NIPC),  [B  -  K  / 
(c°t  -  Lfi^)]  must  be  negative. 

B.  The  c^  elastic  anomaly. 

Within  the  framework  of  the  Landau  theory,  an  effective  elastic 


constant  is  calculated  by  a  set  of  equations 
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3G/3n  *  0;  3G/3e1  -  X1  ■  c^e^, 


(13) 


11 


Using  the  prescription  of  eq.  ( 1 3 )  and  eq.  (8)  gives 


3G/3n  -  (A(T)  -  Ke1  -  Le^Jn  +  Bn3  +  Dn5  -  0 


(14) 


and  3G/3ei  =  (c,,  -  Ln2)  e1  -  1/2Kn2  -  X1 . 


(15) 


In  the  high  symmetry  (high  temperature)  phase  the  strain  and  the 
meter  possess  only  fluctuations,  hence  in  eqs.  (14)  and  (15)  n  * 
6e^ .  Employing  the  standard  procedure,  the  c1 1  elastic  constant 


order  para- 
fin  and  e^  * 
becomes 


c 


11 


for  T  >  T 

o 


in  which  is  the  background  elastic  constant  which  linearly  increases  with 
decreasing  temperature. 

In  the  low  symmetry  (low  temperature)  phase  n  and  e1  possess  spontaneous 
values  n  and  given  by  eqs.  ( 1 1 )  and  (12).  The  order  parameter  and  strain 
can  then  be  described  by  small  fluctuations  fin  and  fie1  around  the  equilibrium 
values  and  ,  and  thus 


"  ■  * 4  Sn  >  ei  •  5i  *  *v 

Eqs.  (14)  and  (15)  can  be  written  to  linear  order  in  fluctuations  as 
(A(T)  ♦  5BfS2  ♦  5D?14  -  Ke1  -  Le2)6n  -  (K*  ♦  2L81fl)«e1  -  0 


(16) 


in  which  rj  and  ?1  are  the  equilibrium  values  given  by  eqs.  (11)  and  (12).  The 
above  thermodynamic  treatment  which  includes  linear-quadratic  and  biquadratic 
coupling  terms  in  the  free  energy  (eq.  8)  gives  only  discontinuous  downward 
displacement  of  the  elastic  constant  below  the  transition.  Thus,  higher  order 
interactions  do  not  explain  the  experimentally  observed  "softening"  of  c^ 
shown  in  fig.  2.  These  enharmonic  terms  further  seem  to  be  of  little  importance 
even  close  to  the  transition  where  critical  fluctuations  in  n  become  large. 

These  calculations  indicate  that  the  explanation  of  the  critical  behavior  of  the 
c1 1  elastic  constant  should  be  sought  elsewhere;  in  particular,  the  dynamics  of 
the  system  must  be  considered. 

C.  Other  mechanisns. 

As  mentioned  above,  the  linear-quadratic  and  biquadratic  couplings 
between  the  order  parameter  and  the  zone  center  elastic  modes  do  not  affect 
the  elastic  properties  within  the  Landau  theory  in  the  high  temperature  phase. 
The  situation  in  the  crystal  is  more  complicated  for  the  case  of  NIPC  and 
other  mechanisms  which  can  at  least  in  part  contribute  to  the  observed  anomaly 


13 


of  the  [lOO]  LA-phonon  should  be  mentioned.  Among  these  are  (l)  enharmonic 
k-dependent  interactions  between  the  [lOO]  LA-mode  and  fluctuations  of  the 
order  parameter,  (2)  a  soft  zone  center  optical  mode  which  couples  to  the 
acoustic  mode,  (3)  a  mean  field  Landau-Khalatnikov  contribution,  and  (4)  dis¬ 
persion  of  the  elastic  constant. 

2 

Anharmonic  coupling  to  critical  fluctuations  of  the  order  parameter  <n  > 

may  lead  to  a  rounding  of  the  discontinuity  in  c^  ^  for  T  >  Tq  and  generate  an 

additional  contribution  in  the  low  temperature  phase. This  mechanism 

accounts  partially  for  the  critical  behavior  in  chloranil1^,  terbium  molybdate^ 

22 

and  ABO^  perovskites.  A  soft  optical  mode  at  the  Brillouin  zone  center  seems 

to  be  one  of  the  possible  explanations  assuming  the  coupling  of  such  a  mode  to 

the  Brillouin  zone  boundary  mode  during  the  transition  exists.  Coupling  between 

strains  and  temperature  dependent  Raman-optical  modes  provides  a  satisfactory 

1 2 

explanation  for  the  elastic  behavior  in  s-triazine  and  the  dependence  of  the 

10  23 

qTA-mode  in  chloranil.  A  mean  field  Landau-Khalatnikov  contribution  asso¬ 
ciated  with  the  dispersion  of  the  elastic  constant  can  contribute  to  the  elastic 
behavior  in  the  low  temperature  phase  as  a  rounding  and  restoring  of  the  elastic 

constant.  This  mechanism  is  known  to  account  partially  for  the  behavior  in 

24 

layered  perovskite  compounds. 

Although  the  above  mentioned  mechanisms  can  all  contribute  to  the 
anomaly  of  the  c^  elastic  constant,  especially  in  the  vicinity  of  the  phase 
transition,  they  cannot  entirely  account  for  the  observed  behavior.  Thus  we 
find  them  possible  but  not  likely.  As  will  be  argued  in  the  next  section,  the 
most  reasonable  explanation  of  the  behavior  of  the  [lOO]  LA-phonon  in  NIPC  is 
based  on  a  dispersion-relaxation  mechanism. 

D.  Dispersion  of  the  elastic  constant  and  dynamical  mechanisms. 

One  possible  explanation  of  the  critical  behavior  of  c^  is  that  the 
elastic  constant  exhibits  a  frequency  dependent  relaxation  which  undergoes 


critical  slowing  down  in  the  vacinity  of  the  phase  transition  (ca.  T  *  137 
40  K).  We  propose  a  semiquantitative  interpretation  of  the  phase  transition 
in  NIPC  based  on  this  assumption.  Since  no  ultrasonic  or  neutron  scattering 
data  are  available  for  NIPC,  the  detailed  mechanistic  dynamics  cannot  be  treated 
in  much  detail.  Measurements  of  the  [lOO]  LA-mode  velocity  by  a  piezoelastic 
"resonance-antiresonance"  method^’ ^  give  a  value  of  2150  ms  ^  at  room  temp¬ 
erature.  The  velocity  is  almost  30%  lower  than  the  velocity  obtained  by  Bril- 
louin  scattering  (2953  ms-^).  In  the  former  method  the  velocity  is  measured  at 
the  frequency  of  a  standing  acoustic  wave  (-90  kHz)  whereas  in  Brillouin  scatter¬ 
ing  the  velocity  is  measured  at  a  frequency  of  -10  GHz.  The  difference  in  the 
sound  velocities  measured  by  these  methods  is  definitely  outside  the  contribu¬ 
tion  from  usual  dispersion  effects,  typically  less  than  3%,  and  certainly  much 
larger  than  the  experimental  errors.  Moreover,  the  temperature  dependence  of 
the  sound  velocity  obtained  by  both  methods  given  in  fig.  8  shows  a  slight 
suppression  of  the  critical  behavior  in  Brillouin  scattering.  In  order  to  con¬ 
firm  the  difference  in  the  sound  velocities  at  room  temperature,  we  have  addi¬ 
tionally  performed  ultrasonic  pulse-echo  experiments  at  a  frequency  of  10  MHz. 

The  velocity  of  the  [lOO]  LA-mode  at  10  MHz  is  2250  ms”^ . 

The  above  experimental  findings  clearly  suggest  that  dispersion/ 
relaxation  and  its  critical  behavior  at  the  phase  transition  are  important  for 
the  elastic  properties  of  NIPC.  Moreover,  the  change  of  the  attenuation  at  the 
transition  indicates  that  some  dynamical  relaxation  process  is  involved  in  this 
behavior.  One  possible  mechanism  for  such  a  relaxation  effect  would  be  a  dipolar 
orientational-relaxation  process.  From  the  standpoint  of  a  "ferroelectric" 
phase  transition,  the  anomaly  in  NIPC  could  be  represented  as  an  order-disorder 
transition  of  the  permanent  dipole  moments  of  the  individual  molecules.  This 
mechanism  remains  in  agreement  with  the  gradual  freezing  of  the  librational 
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motions  of  molecules  observed  in  the  x-ray  measurements  and  with  the  temper- 


4 

ature  dependence  of  the  spontaneous  polarization. 

The  most  straight-forward  way  to  treat  the  frequency  dependence  in  the 
high  temperature  phase  is  to  assume  a  Debye  relaxator  form  for  the  elastic 
constant : 


(20) 


A  good  fit  of  the  experimental  results  in  fig.  2  to  this  equation  may  be 

Q  O 

obtained  with  parameters:  c1 1  (®)  ■  15.11  x  107  Nm  ,  Ac1 1  ■  0.61  x 
109  Nm-2,  Tq  ■  135.5  K,  p  *  1.01;  v  ■  0.44,  (wt  *  1  and  a  =  1 .5  x 


In  the  low  temperature  phase  a  restoring  of  the  c^  elastic  constant 


is  expected  as  temperature  decreases  from  T  .  The  increase  in  c.,  can  be 
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expressed  in  a  manner  similar  to  the  critical  softening  of  eq.  (20) 


(21) 


-v 

T  -T 

•  '  o 

in  which  t  -  t  ( — ^ )  .  The  parameters  which  give  a  good  fit  of 

°  o  , 

eq.  (21)  to  the  observed  behavior  in  the  low  temperature  phase  are  c^(®)  * 

17.19  x  109  Nm-2,  Ac'  -  3-96  x  10-5  Nm"1 ,  T  -  140.8  K,  o'  -  0.43,  (on)2  * 

I  1  o  o 

1  and  v  *  0.27.  If  one  expands  eq.  (20)  and  (21)  in  a  power  series  in  (T-7q)/ 
Tq  and  (Tq-T)/To,  respectively,  eq.  (4),  previously  used  for  fitting  the  aata, 
obtains.  The  fits  of  eqs.  (20)  and  (21)  to  the  data  are  not  completely  unique 
because  of  the  number  of  available  parameters.  Consequently,  a,  Tq,  (®) 
and  c  ,  ( “)  are  chosen  to  correspond  to  the  appropriate  values  obtained  by 
fitting  the  phenomenological  eq.  (4)  to  the  data.  Moreover,  since  for  Brillouin 


scattering  typically  tqio  -  1  near  the  transition,  t^uj  £  1  in  the  above  procedure 
Nonetheless,  the  values  of  Ac^ ,  pand  v  are  not  completely  determined  and 
siminlar  quality  fits  can  be  obtained  for  small  (_♦  20%)  variations  of  all  three 
parameter  values. 

Assuming  a  Debye  form  for  the  dispersion,  the  relaxation  time  may  be 
obtained,  as  shown  in  the  case  of  barium  sodium  mobate,  by  combining  the 

temperature  dependence  of  c^  from  fig.  2  with  the  increase  of  attenuation 
displayed  in  fig  6.  Employing  the  theory  of  a  complex  elasticity  with  a  single 
relaxation  time,^  the  elastic  constant  may  be  written  as 

C1 1  (“ )  *  “  1  -  0)T  (T)  *  (22) 

Separating  eq.  (22)  into  real  and  imaginary  parts  yields 
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By  combining  these  two  equations  one  obtains  a  very  simple  expression  for  the 
relaxation  time 


t(T) 


ril<T>-r, 
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2ttP 

2  * 
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(25) 


Equation  (25)  is  valid  even  if  the  relaxation  time  and  coupling  are  wavevector 
dependent,  since  it  involves  data  from  a  single  frequency  u>  and  wavevector  q . 
fig.  9  shows  the  inverse  of  the  relaxation  time  as  a  function  of  temperature. 


The  line  in  the  high  temperature  phase  is  of  the  form 


T_1(T,q) 


-1  ,  -1  ,  . 
t1  +  t2  (q)  . 


(26) 
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The  fit  to  this  equation  gives  t1  =  1.16  x  10  s  and  ^ 
the  low  temperature  phase 


1 .08  x  10“10s.  In 


(27) 


•  _i  ?  '  -10 

with  values  of  the  coefficients  *  2.4  x  10  s  and  *=  1.07  x  10  •  The 

above  temperature  dependences  of  the  relaxation  times  satisfy  the  prediction 
of  a  mean  field  critical  slowing  of  the  order  parameter  fluctuations  in  the 
vacinity  of  Tq. 

V.  Conclusions. 

The  elastic  properties  of  NIPC  in  the  temperature  range  90  K  -  295  K 
have  been  investigated  by  Brillouin  scattering.  In  particular,  the  nonferroic, 
first  order  phase  transition  at  ca.  137  K  has  been  examined.  Our  findings  can 
be  summarized  as  follows: 

1.  The  c^  elastic  constant  associated  with  the  [ 1 00 ]  LA-phonon 
exhibits  an  anomaly  extending  ca.  40  K  above  and  below  the  phase  transition  at 
137  K.  The  anomaly  has  the  form  of  a  softening  of  the  c^  elastic  constant; 

falls  to  a  minimum  value  in  the  transition  region.  The  halfwidth  at  half 
maximum  of  the  [lOO]  LA-phonon  also  possesses  anomalous  broadening  around  the 
transition  temperature. 


The  0^2*  c^f  c ^  and  most  likely  also  c^  elastic  constants 


exhibit  a  step-like  anomaly  at  ca.  137  K  characteristic  of  a  first  order  transi¬ 


tion.  Practically  no  anomalous  behavior  is  observed  for  the  c ^  elastic  constant. 
Only  weak  anomalies  for  the  c^i  end  C2J  elastic  constants  are  found. 


3«  Sound  velocity  diagrams  in  the  three  crystallographic  planes  are 
calculated  at  room  temperature  using  the  values  of  the  elastic  constants  deter¬ 
mined  in  our  experiments. 

In  order  to  explain  the  anomalous  behavior  of  the  elastic  constant 
two  possibilities  are  considered:  1)  an  elaboration  of  the  Landau  theory  with 
the  inclusion  of  the  linear-quadratic  and  biquadratic  couplings  between  the  zone 
boundary  order  parameter  and  zone  center  elastic  strain;  and  2)  dispersion  of 
the  elastic  constant  combined  with  a  temperature  dependent  relaxation  process. 

The  c^  anomaly  cannot  be  understood  within  the  framework  of  the  Landau 
theory  which  can  only  account  for  the  downward  discontinuous  change  of  c^  at 
the  transition.  The  anomaly  can  be  understood  by  assuming  dispersion  of  the 
elastic  constant,  possibly  due  to  some  relaxation  process.  This  dynamical 
approach  is  justified  by  experimental  findings:  a  large  difference  in  the 
velocity  of  the  [ 1 00 ]  LA-phonon  measured  by  experimental  techniques  employing 
different  frequency  domains;  and  suppression  of  the  critical  behavior  at  high 
frequencies  (“10  GHz).  In  the  framework  of  this  latter  approach  the  experimen¬ 
tal  data  are  reproduced  using  a  power  law  equation  including  a  term  due  to  the 
assumed  Debye  relaxation  behavior.  The  temperature  dependence  of  the  relaxation 
time  is  calculated  by  combining  the  temperature  dependences  of  the  elastic 
constant  and  attenuation.  This  behavior  is  found  to  satisfy  the  mean  field 
dependence  characteristic  of  critical  slowing  down. 
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TABLE  I 


Solutions  of  the  equation  of  motion  for  NIPC.  Only  modes  observed  in 
the  experiments  are  included  in  the  table.  [See  eqs.  2  and  3  in  text]. 


Phonon 

Direction 

Polarization 

Xi  =  pv2 

L  [100] 

xi  ■  cn 

[100] 

T  [OIO] 

X2  =  C66 

T  [001 ] 

X3  =  C55  +  e15/e11 

L  [OIO] 

X4  ‘  C22 

[oio] 

T  [lOO] 

X5  *  C66 

[001  ] 

L  [OOI] 

Xb  *  C33  +  e33/e  33 

T  [lOO] 

X7  -  C55 

qL  [l10] 

*8  •  lC11  ' 

f  C22  +  2C66  +  ^C11  " 

C22)2  *  1(C12  * 

[110] 

qT  [l10] 

X9  ‘ 

_  II  _ 

_  M  _ 

T  [OOl] 

X10  '  '/ 2  (C44 

+  C55  +  (e15  +  e24)?/( 

£11  +  e22^ 

[101] 

qL  [101] 

X11  ■  V4  |C„ 

*  .  2C„  *  t(c„  - 

C„)2  .  4<C„ 

* 

qT  [101] 

X12  “ 

_  II  _ 

_  n  _ 

[oil] 

qL  [Oil] 

x„  •  1/4  |C2J 

+  C33  *  2C44  +  ^(C22  * 

033)  ♦  4(C23 

* 

qT  [Oil ] 

X14  " 

_  «l  ^  _ 

—  H  • 

FIGURE  CAPTIONS 


Figure  1 


Figure  2 


Figure  3 


Figure  4 


Figure  5 


Polar  plots  of  the  sound  velocity  of  different  acoustic  inodes 
in  the  ab  (a),  ac  (b),  and  be  (c)  planes  of  NIPC,  obtained  using 
the  nine  independent  elastic  constants  listed  in  Table  II.  The 
lines  are  calculated  using  the  general  solutions  of  the  equation 
of  motion,  eq.  (2).  Particular  solutions  in  the  directions  of 
the  crystallographic  axes  and  diagonals  are  given  by  relations 
from  Table  I.  The  units  are  1(r  ms 

Temperature  dependence  of  the  c1  elastic  constant.  (o), 

(x)  correspond  to  two  different  experimental  runs  performed 
using  the  samples  of  different  orientations.  The  two  runs  are 
normalized  at  190  K.  The  lines  represent  fits  to  the  phenomeno¬ 
logical  eq.  (4)  with  the  values  of  the  parameters  given  in  the 
text. 

Temperature  dependences  of  the  c^2  and  c ^  elastic 
constants:  (o)  and  (x)  represent  different  experimental  runs. 

Temperature  dependences  of  the  c^,  c ^  and  c ^  elastic 
constants:  (o)  and  (x)  correspond  to  two  independent  experi¬ 

mental  runs. 

The  temperature  development  of  the  NIPC  Brillouin  spectrum  due 
to  [lOO]  phonons.  The  free  spectral  range  is  21.8  GHz.  The 
longitudial  and  transverse  modes,  labeled  L  and  T,  respectively, 
correspond  to  the  Rayleigh  peak  marked  R.  The  spectra  are  taken 
in  right  angle  scattering  (e  *  90°)  in  triple-pass. 


Figure  6 


Temperature  dependence  of  the  Brillouin  halfwidth  at  halfmaximum 
for  the  Brillouin  feature  corresponding  to  the  [ 1 00]  LA-mode  in 
fig.  5.  The  halfwidth  was  not  deconvoluted  with  respect  to  the 
width  of  the  Rayleigh  peak,  (o)  and  (x)  correspond  to  the  two 
independent  experimental  runs.  The  lower  points  represent  the 
temperature  dependence  of  the  Brillouin  intensity  of  the  same 
Brillouin  feature  taken  as  the  area  underneath  the  peak. 

Figure  7  Temperature  dependence  of  the  effective  elastic  constants 

corresponding  to  [ 1 1 0]  (a),  [ 1 01 ]  (b)  and  [oi 1 ]  (c)  acoustic 
phonons  specified  in  Table  I.  (o)  and  (x)  correspond  to  differ¬ 
ent  experimental  runs. 

Figure  8  Temperature  dependences  of  the  sound  velocity  of  the  [ 1 00 ]  LA- 

mode  obtained  by  Brillouin  scattering  (crosses)  at  a  frequency 
12  GHz  and  "resonance-antiresonance"  method  (circles)  at  a  frequency 
90  KHz.  The  scales  are  slightly  shifted  with  respect  to  one 
another  for  easy  comparison. 

Figure  9  Temperature  dependence  of  the  inverse  relaxation  time  associated 

with  the  elastic  constant.  Solid  lines  are  the  least 
squares  fits  to  eqs.  (24)  and  (26)  in  the  high  and  low  tempera¬ 
ture  phase,  respectively. 
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